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1. Introduction 

In this paper we will examine the thin string limit of solutions of the field equations for 
an infinite length straight Cosmic String [described by a complex scalar field coupled 
to a U(l) gauge field] coupled to gravity. The thin string limit has been examined by a 
number of authors [1-4] and it is generally agreed that in the thin string limit the metric 
describing an infinite length static Cosmic String becomes that of a conical spacetime 

ds^ = dt" - dr^ - A\^d(f)^ - dz" (1) 

where (f) is the usual 27r-periodic coordinate of cylindrical polar coordinates. If one 
changes to a new angular variable = Acj), where (f) e [0, 27rA], the metric becomes that 
of Minkowski space so we see that (1) represents a conical metric with angular deficit 
A0 = 27r(l — A). If one also makes the weak-field approximation then one can also show 
(see e.g. Vilenkin [|l|) that the angular deficit A0 is related to the mass per unit length 
/i of the the string in a simple way: 

A0 = Svr/i (2) 

(where we are using units in which h = c = G = 1). 

This result is consistent with the study of "simple line sources" in general relativity 
made by Israel [Q. He showed that if one has a 1-parameter family of spacetimes with 
the matter tending to a simple line source in the limit then, providing the radial stress 
becomes negligible compared to the mass density, the energy-momentum tensor becomes 
that of an infinitely thin wire with distributional stress energy. This analysis was applied 
by Linet [^ to a critically coupled Cosmic String (i.e. a string for which the scalar and 
vector fields have equal masses) who showed that equation (2) is satisfied in this case. 
However as pointed out by Garfinkle [^, in the non-critical case Israel's condition is 
violated and there fails to be a simple relationship between A0 and in the thin string 
limit. Futamase and Garfinkle [Q argue that because the relationship between A0 
and /i depends upon the nature of the matter which makes up the string there is in 
general no well defined thin string limit for Cosmic Strings in general relativity. This 
interpretation agrees with the analysis of conical spacetimes by Geroch and Traschen 0], 
who showed that different 1-parameter families of spacetime with a conical spacetime 
as the limit could have different limiting mass per unit length. Note however that if one 
puts certain additional conditions on the 1-parameter families of spacetime then one 
does recover formula (2) (see for example Wilson |Q and Clarke et al [Q). 

The difficulty of having a well defined thin string limit in general relativity comes from 
the non-linear nature of the field equations. For a linear theory such as electromagnetism 
a consistent framework for describing concentrated sources is provided by Schwartz's 
theory of distributions. This is for two reasons. Firstly the equations are linear with 
respect to sources and fields so the field equations are well defined for distributional 
sources. Secondly for smooth sources close to a point charge (in the topology on the 
space of distributions) the field produced is close to that produced by a point charge. 
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It is this second feature which means that in electromagnetism the detailed internal 
structure of a concentrated source is unimportant and only the distributional charge 
density that it defines in the limit is significant. 

In this paper we will consider the thin string limit of Cosmic Strings by making 
use of generalised functions belonging to the simplified Colombeau algebra Qs- The 
construction of the simplified Colombeau algebra is outlined in section 2 but the key 
point to note here is that it provides a consistent setting which allows one to multiply 
distributions in a non-linear theory such as general relativity. Furthermore if two 
generalised sources are close (in the appropriate topology) then the corresponding 
metrics (satisfying the appropriate boundary conditions) are also close. Although 
Qs contains distributions D'(M^) as a vector subspace, the algebra has a much richer 
structure because it is able to encode non-linear information. Thus one can use these 
algebras to model the thin string limit of a Cosmic String coupled to gravity because 
one is able to retain certain features of the source in the limit. This explains why one 
can use them to define the thin string limit even in case when there is not a simple 
relationship between A0 and fi. 

The plan of the paper is as follows. In section 2 we outline the construction of the 
simplified Colombeau algebra. In section 3 we illustrate the ideas by considering the 
thin string limit of a static Cosmic String in Minkowski space. In section 4 we consider 
the thin string limit of a static Cosmic String coupled to the gravitational field and in 
section 5 examine the relationship between the angular deficit and mass per unit length. 
In section 6 we introduce the conical approximation and in section 7 we briefly consider 
the case of a dynamic Cosmic String coupled to gravity. 



2. A brief review of Colombeau theory 

In this section we give a brief introduction to the Colombeau algebras. For more details 



see 



Ton and [0. The basic idea is to consider generalised functions as 1-parameter 
families of smooth functions {/e}. Our basic space will thus be 

f (M") = {{/,} : < e < 1, /, G C°°(M")}. (3) 

We may represent a smooth function by the constant family fe{x) = f{x), but for a 
function of finite differentiability we may obtain a family of smooth functions by taking 
the convolution with a suitable smoothing kernel or molifier $ 

U^) = \ I f{ym{y-x)/e)<ry. (4) 

However for the case of a smooth function we would like both representations to be 
equivalent. Using a Taylor series expansion to compare the difference between these 
expressions we are lead to define two representations to be equivalent if they differ by 
a 'negligible function' which is defined as a 1-parameter family of functions which on 
any compact set vanishes faster than any given positive power of e. Since we are trying 
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to construct a differential algebra wc also require the derivatives of / to satisfy this 
property and the resulting set J\f to be an ideal. Unfortunately A/" is not an ideal in 
£^(]R"), but by restricting the space to 'moderate functions' Sm{^^'') which grow no faster 
than some inverse power of e, one does have an ideal and we may define the differential 
algebra Qs as the quotient. 

Definition (Moderate functions) 

£:m(M"H{/6} : CC W, Va e N^, 3p e N, 3?7 > 0, 3C> such that 
Sup|r'"/,(x)| ^ Ce~P for < e < 77} . 



Note K CC indicates that K is compact and we have also employed the standard 
multi-index notation for D"/. 

Definition (Negligible functions) 

J\f{Wy{{f,} : yK CC M", Va G N^J, Vg e n,3r] > 0,3C > such that 
Sup ID"/, (a;) I ^ Ce'^ for < e < 77} . 



Definition (Simphfied algebra) 

^,(R") = £:m(m")/a^(k"). 

Thus a non-linear generalised function is represented by a moderate sequence of smooth 
functions modulo a negligible sequence. The space Sm(^") is a differential algebra with 
pointwise operations and since the space of negligible functions is a differential ideal, Q 
is also a commutative differential algebra. As remarked earlier one of the advantages 
of the Colombeau approach is that one may frequently interpret the results in terms of 
distributions using the concept of association or weak equivalence. An element [/,] of Q 
is said to be associated to a distribution T e if a representative {/g} in Sm satisfies 

V(/)eP, lim/" f,{x)(f){x)(rx^ {T,(f)) (5) 

and we then write [/g] fti T. Note that not all elements of Q are associated to 
distributions. 

More generally we say [f^] ~ [g^] if 

V(/) e V, lim [ f,{x)(f){x)(rx = lim [ ^g(x)0(x)d"x. (6) 

Association is an equivalence relation which respects addition and differentiation. It also 
respects multiplication by smooth functions but by the Schwartz impossibility results 
cannot respect multiplication in general. 

The algebra presented above is the simplest of the Colombeau algebras. Furthermore 
the notion of smooth function and the growth conditions in the definitions of moderate 
and negligible functions only depend upon the differential structure. This allows the 
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simplified algebra to be readily generalised to arbitrary manifolds. However it does 
suffer from the disadvantage that the embedding of functions and distributions is 
not canonical. Thus one has to appeal to mathematical or physical arguments outside 
the theory to justify a particular representation of a non-smooth function. For the case 
of the field equations for a Cosmic String coupled to gravity, we will see in the next 
section that there is a natural scaling of the coupling constants which can be used as 
the Colombeau parameter. In other cases one can embed by taking the convolution with 
an element of some special function space and only use general properties of elements 
of the space in proving results. However this is not always satisfactory and in certain 
situations the full Colombeau algebra is preferable. This involves enlarging the algebra 
substantially by making the generalised functions depend upon the molifier. This allows 
one to construct a canonical embedding of distributions but at the expense of problems 
with coordinate invariance. A manifestly coordinate independent construction of the 



full algebra is given in [12 



3. The thin string limit in Minkowski space 

We consider Cosmic Strings which are described by a complex scalar field $ which is 
coupled to a vector field A^. The equations of motion for these fields may be derived 
from the Lagrangian 

L = D^^W^-V{\^\) - ^F^uF^" (7) 

where = = ieAfj^ is the gauge covariant derivative with coupling constant e and 
Ffj.u = V^v4y — V^Af^. We write $ in the form $ = -^Se^"^ and consider the case of a 
potential given by \^ = 4A(S'^ — r/^)^. Then the Lagrangian becomes 

Here A is a coupling constant and rj^ is twice the vacuum expectation value of the 
scalar field. For the above Lagrangian the mass term for the scalar field is given by 
m| = 8A?7^ while that for the vector field is given by mp = e^rj^ so that the mass terms 
are equal and the string is critically coupled when = 8A. The energy-momentum 
tensor corresponding to the Lagrangian @ is given by 

T^u = V^SV.S + S^V^tfj + eA^){V.^ + eA,) - F/F,r - g^^uL. (9) 

We now look for solutions describing an infinite length static cylindrically symmetric 
Cosmic String in Minkowski space. Following Garfinkle we therefore assume that 
the scalar and gauge fields have the form 

$ = S{p)e''"^ (10) 

T) 

/l, = -[P(p)-l]V,0 (11) 

where (t, z, p, 0) are the usual cylindrical polar coordinates for Minkowski space. The 
number n is the winding number for the string which for simplicity we will take to be 
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one in the remainder of the paper. With this choice the equations of motion are given 

by 



d / dS\ 



^ dp VdpJ 



dp\ dp / 



(12) 
(13) 



The requirement that the solution has finite energy imphes that 

S{p) rj, as p CO (14) 

P(p) ^0, as p ^ CO. (15) 

On the other hand at the centre the solution must smoothly attain the symmetric state 
which requires 

S{0) = (16) 
P(0) = 1. (17) 
The fields describing the Cosmic String are therefore solutions of (|T2|), (^) satisfying 

With the fields taking the form (|10|), (|1T]) the energy-momentum tensor is diagonal and 
is given in terms of S and P by 



where 



Pp 









'dsy 


-Pz = 5 


[( 


. dp J 


1 


-^dS\ 


2 




2 






p2 



:2 p2 



1 /C/P\2 



e^p^ \dp ) 



-(— I + — ^ 2X{S^-ri^f + ( — 1 



\dp ' 



P 



e^p^ \dp / - 



The mass per unit length in Minkowski space is therefore given by 

crd^x = 27r 



p 



R2 



crpdp 



(18) 

(19) 
(20) 
(21) 

(22) 



where a is given by (|TUp. 

Although the fields extend to infinity, they decay exponentially fast. For large p the 
asymptotic behaviour is given by 



P(p) ~ pKi{er]p) ~ a / ^P"'e 



7r 



l/2^-e,?p 



S{p) T]- K(i{V8Xrip) T]- 



TT 



-P 



-1/2-VSXvp 



(23) 
(24) 



V32Xr] 

so that cr also decays exponentially fast. We therefore define an "effective radius" p for 
the Cosmic String by requiring that 90% of the total mass per unit length is within this 
radius and hence 



|.p g i^oo 

(^{p)p^p = Yo io 



(25) 
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Note that the effective radius is roughly the Compton wavelength of the scalar field 
which for typical values of the coupling constants is of the order of 300 Planck lengths. 
If we define r = \^rip, then numerical investigations show that r is very close to 1 as 
e, A and r] vary over several orders of magnitude. 

Thus in practice most of the matter is highly concentrated on the axis of the string. For 
this reason it is natural to try and model a Cosmic String by taking a thin string limit 
in which all the matter is concentrated on the axis. However as noted by a number of 
authors (see eg 0) there are problems in simply using a 5- function to model the source. 
The approach adopted in this paper is to allow e A and rj to depend upon a parameter 
e so that 

e = e(e), A = A(e), V = v{^), 

with the physical values given by e = 1 and the thin string limit given by letting e tend 
to zero. We then consider the corresponding 1-parameter family of solutions >S'e(p) and 
Pe{p) of (|r2p-(|T^) which we regard as elements of £m, whose equivalence classes [S^] and 
[Pg] are elements of Qg. Since Qg is a differential algebra we may compute the energy- 
momentum tensor (and other physically interesting quantities) within Qg. Whether 
these quantities have an interpretation in terms of the usual theory of distributions will 
depend upon whether they are associated to a conventional distribution. 

We now consider the properties we require of the thin string limit in order that it 
retains the important features of the the physics. In order that e represents the 
thin string limit we require that the effective radius p(e) defined by (|25D scales linearly 
with e, while the mass per unit length /i(e) is independent of e. We also require that 
the relative strength of the coupling between the scalar and vector field given by e^(e) 
compared to the self-coupling of the scalar field given by A(e) is independent of e. We 
therefore try and choose the scaling for e(e), A(e) and ri{e) so that 

Pie) = epil) = ep (26) 

/i(e)=/i(l)=/i (27) 
e\e)/\{e)=e\l)/\{l) = a. (28) 

In fact conditions (p6|) -(p8l) uniquely determine the e dependence of e(e), A(e) and 77(e) 
to be given by 

e(e) = e/e (29) 

A(e) = A/e^ (30) 

ri{e)=r]. (31) 

To see this we start by observing that after a suitable rescaling the only significant 
parameter in (^2[) and (|13[) is the ratio of the coupling given by a. Thus a solution to 
(0) and ([ID satisfying (|13)-(|13) is given by 



Sip) = vie)Xiy^vie)p) (32) 
P{p) = Yi^vie)p) (33) 
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where X{r) and Y{r) are solutions of 

d / ^dX\ ^ x[A\r\X^ - 1) + Y'] (34) 

dr\ dr / 
satisfying the boundary conditions 

X(0) = 0, hm X(r) = 1 (36) 

r(0) = 1, hm Y{r) = 0. (37) 

r— >oo 

Substituting for S{p) and -P(p) using (^) and ( ^3]) in the expression ([l9|) for a we find 
that 



^e(p) = A(e)r/^(e)Mv^r^(6)p) (38) 

where 

h{r) = d[X'{r)r + \x'{r)Y'{r) + 2[X\r) - 1]' + —.[Y'ir)]'}. (39) 
The corresponding mass per unit length is then given by 

POO 

/i, = 277 / a^{p)pdp 
Jo 

/"OO 

= 271 r]'^{e)h(r)rdr 
Jo 

= rf{e)M (40) 

where 

/■OO 

M = 2n h{r)rdr. (41) 

Since M is independent of e we see that for the mass per unit length to be independent 
of e we also require that ri{e) is independent of e so that 

vie) = r/(l) = r]. (42) 

We next introduce a scaling factor s(e) defined by 

. £M = £(£) (43) 

e(l) e 

Then in order to make e^(e)/A(e) independent of e we see that A(e) must scale according 
to 

A(e) = s\e)X. (44) 

We now calculate how the effective radius depends upon s(e). We first let p{p) be the 
total mass per unit length within a radius p so that 

f 

fJ,e{p) =271 (Te(x)xdx 



2717]"^ / h(r)rdr 



H{s{e)p) (45) 
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where H{r) is given by 

H{r) = j h{s)sds. (46) 

Now when p = we require 

9 9 
We) = Y^f^e = — 

so that we want 

o 9 
H{s{e)p,) = Y^/^ = const. 

But by requirement (^Bf) p^ = ep so that we need 

H{s{e)ep) = const. (47) 

which imphes that s(e)e is constant (since H is increasing) and since we require s(l) = 1 
we must have 

s{e) = 1. (48) 
We have therefore estabhshed the following proposition. 

Proposition 1 The e- dependence of e{e), A(e) and 77(e) such that a solution of the 
field equations for a Cosmic String in Minkowski space satisfies conditions ^dj)-^^) is 
uniquely given by 

e{e) = e/e, A(6) = X/e\ r/(e) = 

and the corresponding solutions of ^1% ) and ( ff^ j satisfying are given by 

Se{p) = vXiVXvp/e) (49) 
P,(p) = r(v^r/p/e) (50) 

where X{r) and Y{r) are solutions of ij^2) and iWW satisfying the boundary conditions 



We next calculate the corresponding value of the energy density <7e{p)- Substituting 
(^) and ( |50D into the expression for a (|T9|) we find 

o.(p) = ^<^) = >(pA). (51) 

We now show that <Je{p) represents a 5-function within Qs- Because of the coordinate 
singularity at the origin in polar coordinates we first transform to Cartesian coordinates 
and define ae{x,y) = ae{^/x^~-\-y^) and de{x,y) = {1/ p)ae{x,y). Then 

/ de{x,y)dxdy = - ae{x,y)dxdy 



— [ (^eip)pdp 

P Jo 



Pe 
P 

1 (since Pe = p Ve). (52) 
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Also 



A'' 



= —di{x/e,y/e). 
Thus ^^(x) is a 1-parameter family of smooth functions satisfying 



(53) 



(i) 



. 1 



d 



Using the results of section 2 we see that [de(x)] represents a 5-function in Qg- Note that 
although the concept of delta-function in Qs is not unique, the various delta-functions 
are all associated to the usual distributional 5-function so that 

[a,ix,y)]^fiS^'\x,y). (54) 

Using a similar argument one can formally show that the principal pressures Pp and 
are proportional to 5-functions 

[p;{x,y)]^a6^'\x,y) (55) 
[plix,y)]^b5^'\x,y) (56) 



where 



27T 



2tx 



Ppip)pdp 

) 

p^{p)pdp. 



(57) 
(58) 



However because of the degeneracy of the polar coordinates at the origin the meaning 
of (^) and (|56|) is unclear. We therefore transform to Minkowski coordinates {t, x,y, z) 
and find that the corresponding distributional energy-momentum tensor T^^, is given in 
these coordinates by 





T, 



/ 










\ 







c 
















c 





\ 











p J 



6^^\x,y) + d 



/O 






2xy 







"ixy 

9 9 

y^—x^ 
x'^+y^ 





S^'\x,y) 



0\ 





0/ 

= COS 29 and 



(59) 



^xy 
x'^+y'^ 



sin 26 are 



where c = |(a + 6) and d = ^{c + d). We note that ^■^_^y-^ 
bounded but have a directional dependent limit. Indeed the functions cos 2^ and sin 2^ 
are precisely the directional dependent functions that appear in the conical metric when 
written in Cartesian coordinates (see for example |^). More importantly in the classical 
theory of distributions one can only multiply elements of P'(]R^) by smooth functions 
so that the entries in the second matrix in ( ^9]) do not define elements of P'(R^). It is 
not hard to extend the notion of multiplication in X>'(R^) in some ad hoc way to allow 
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these sort of products (see for example Oberguggenberger |T^) but instead we prefer 
to calculate the energy- momentum tensor within Qg. To do this we first calculate the 
(e-dependent) energy-momentum tensor in Minkowski coordinates for the 1-parameter 
family of solutions given by (^) and (|50D. This tensor is regular (even on the axis) 
and the components define elements of Sm whose equivalence classes define elements of 
Qg. Thus [Tl^i^{t, X, y, z)] for /x, i/ = 0, ... 3 is a well defined element of Qg. One can then 
show, for example that 

2dxy 
'x^ + y'^ 

where the product on the right hand side is within Qg (and hence well defined). 



(60) 



3.1. Critical coupling 

An important special case of the field equations occurs for critical coupling, where the 
masses of the scalar and vector fields are equal. With our choice of numerical factors in 
the definition of the fields and coupling constants, critical coupling occurs when = 8A 
(or a = 8). In this case one can show that the field equations reduce to the Bogomol'nyi 
equations which in cylindrical symmetry are given by 
dX XY 

— = 61 

ar r 

dY 

= _4r(i_x2). (62) 

dr 



Substitution into ( pil| ) and (^Tj) shows that for the case of critical coupling 

Pp = (63) 
= (64) 

and thus a, 6, and hence c and d vanish. Furthermore in the case of critical coupling 
one may integrate ([191) by parts and use (|6TD and (^) to show that for winding number 
n the mass per unit length is given by 

/i = Tmrf (65) 

which leads to 

[T;,(t,x,i/,^)] ^/i5(2)(a;,Z/)diag(l, 0,0,1). (66) 

Thus in the case of critical coupling the distributional energy-momentum tensor takes 
the form that is usually used in the thin string limit. As we have seen above this is a valid 
choice for a critically coupled Cosmic String where pp and p^ vanish and for such strings 
it is reasonable to describe the energy-momentum tensor using conventional 5-functions 
which do not model any internal structure. However for non-critical coupling the 
energy-momentum tensor is more complicated and cannot be modeled adequately using 
conventional distributions but can be described using non-linear functions belonging to 
the simplified Colombeau algebra. Because of the 'micro-local' structure of elements 
of this algebra one can take the thin string limit but still retain some features of the 
internal structure of the source in the limiting energy-momentum tensor. 
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4. Static Cosmic Strings coupled to gravity 



In this section we examine the physically more important case of a Cosmic String coupled 
to gravity. We again consider the case of an infinite length static Cosmic String, so the 
Lagrangian is still given by (P), with $ and taking the form given by ( ]T0| ) and 



n]), but with the indices in the Lagrangian raised using a general static cylindrically 



symmetric metric rather than that of Minkowski space used in the previous section. 
Such a metric may be written in the form 



ds' = e^'^'^dr' - dp' 



(67) 



This is a special case of the generalised Jordan, Ehlers, Kundt and Kompaneets metric 
(with the twist uj = 0) used by Sjodin et al [0 in analysing dynamic Cosmic Strings. 
The two metrics are related by taking a = 2(7 — ip), b = —2ip and c = 2{ip + /i). Note 
also that a, b and c are slightly different from the form chosen by Garfinkle ||^ in his 
analysis. 

We now define four orthonormal covector fields according to 



T,, 



The energy-momentum tensor is given by 
form of metric (^^ is given by 



4>,i = pe' 



(68) 



which with the choice (p!OD, ([TT|) and the 



where 



-Pz 



g2p2g6 \ dp J . 



dp J 



2X{S' - 
- 2X{S^ 



/dpy 



g2p2g6 \dp J . 



(69) 

(70) 
(71) 
(72) 



\dp J ' p2gb ^ I / • g2p2gf) \ dp J 

Before considering the full field equations we consider a combination of Einstein 
equations which allows us to simplify the metric. We find 

i?;..(r^r- + ^^) = 87r(T^, - i(7^.)(f^r- + PP^) = (73) 
since a = —pz- If we calculate the left hand side in terms of a, b and c we find 



d^ , 
di^^' 

which implies that 
d , 



. d . 



p 2 dp 



0. 



E- 



g-5(^+'=) 



(74) 



(75) 



dp ' p 

where ii^ is a constant. Regularity on the 2;-axis implies that = so that the relation 
([f5|) implies that a = c + const. Without loss of generality we may rescale p and z so 
that a(0) = c(0) = and hence 

a = c Vp. (76) 
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Note that c does not appear in equations ([70|)-(|72D, so the expression for the energy- 
momentum tensor is unchanged by this simphfication to the metric. 

The boundary conditions for the matter fields are the same as those in Minkowski space 
and are given by (p!^-(p^. For the metric, elementary flatness on the axis together with 
the choice a(0) = requires that 6(0) = 0, while the stronger condition of C^-regularity 
implies that a'(0) = and 6'(0) = [0. 



4.1. The thin string limit 

As in the case of a Cosmic String in Minkowski space one may rescale the equations 
to obtain differential equations involving just the ratio a = e^/A rather than e and A 
separately. We therefore introduce the rescaled field and metric variables X, Y, A and 
B, which are defined by 

Sip) = vXiVXvp), (77) 

P(p) = YiVXT^p), (78) 

a{p) = c{p) = A{v^r]p), (79) 

b{p) = B{VXr]p). (80) 

We also follow Garfinkle |Q and further simplify the equations by replacing B{r) by the 
new variable K{r) where 

K{r) = re5(^+^) (81) 

Then X{r), Y{r), A{r) and K{r) are solutions of 

K^(k^) = -87rr^Vf2ir2(X2 - 1)^ --(^V] (82) 

.d^K 



= -Mrfe^[2K\X^ - \f + e^X^F^i (83) 

dr^ 

k4- ( K^) = e^X[AK\X^ - 1) + e^Y^] (84) 
dr \ dr / 

K /e^ dY\ _ ^^Aj^2y (35) 
dr\K dr J 
satisfying the boundary conditions 

X(0) = 0, lim X{r) = 1, (86) 

r(0) = 1, lim F(r) = 0, (87) 

r— >oo 

A{0) = 0, A'{0) = 0, (88) 

K{0) = 0, K'{0) = 1. (89) 

One may now consider the thin string limit by allowing e, A and r] to depend upon e as 
in the previous section. For a string in a curved spacetime we define the effective radius 
p to be the proper radial distance which contains 90% of the mass of the string. Then 
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using similar arguments to the flat space case one can establish the following proposition 



Proposition 2 The e- dependence of e{e), A(e) and 77(e) such that a solution of the 
field equations for a Cosmic String coupled to gravity satisfies conditions ^2di)-^Bi) is 
uniquely given by 

e{e) = e/e, A(e) = A/e^, r^{e) = ^ 
and the corresponding solutions of the field equations satisfying are given by 

Sip) = r/X(v^r/p/e) (90) 
P,(p) = Yiv^r^p/e) (91) 
a,(p) = A{VXr]p/e) (92) 
6,(p) = B{VXr^p/e) (93) 
q(p) = C{VXvp/e) (94) 

(95) 

where r'^e^^^'' = K'^{r)e~^^^^ and X{r), Y{r), A{r) and K{r) are solutions of (jg^j-ffg^j 
satisfying the boundary conditions (^^-((5^. 

The corresponding energy-momentum tensor is given in the (r, p, 0, z) coordinates by 

t;. = ^.(p)^^^ +Pp(p)p;.p: 

where 



and 



A{VXnp/e)':0 C€ _ A{VXrip/e)':l 



(96) 
(97) 



with 



Xr]^^fVXr]p\ 
a.ip) = —a[^), 

f Vxvp \ 



)' 



+ 



Pr 



P4> = 



\dr J ' 



+ 2(X2- 1) 



1 /rfy\2n 



a 



\dr J 
fdX 

\ dr J 



K^\dr J 
dY\^- 



-e"^( — 1 + — — 2{X^-lf + 1 1 



(98) 
(99) 

(100) 
(101) 
(102) 



' aK^KdrJ 

Note that if we transform to rescaled coordinates (t, r, 0, z) with t = y/Xrjr and r = \f\r\p 
[where A = A(l) and ri = ri{l) are the physical values] then the expression for the energy- 
momentum tensor becomes 



= ^e{r)tltl + Wr{r)7lrl + p^,(r)« + pl{rWl 



(103) 



where 

_1_ 1_ _1_ 1_ 

^^.l'^) = — '^('^/e), ^rir) = —Pp{r/e), p'^{r) = —p^{r/e), pl{r) = —p,{r/e), (104) 



The thin string limit of Cosmic Strings 



15 



and 



(PI = Vxv^- 



The above expression may then be used to obtain both the physical values of the energy- 
momentum tensor (by taking e = 1) and the thin string limit (by taking e — 0). 

Before taking the thin string limit we give an expression for the mass per unit length 
(taking e = 1) in terms of rj'^. We start by noting that 



(105) 



H = 2n a{r)K{r)dr 



2v2 



i\2 



Kdr. 



We now complete the square on the (X')2-terms and integrate by parts to obtain 



(106) 



fl = TVI] 



TXT] 



X' 



X' - 



ne^XY\ 2nX'XY n- 



nrj 
nY' 



X' - 



K 

ne^XY 
K 

ne^XY 



+ 



K 

nX'^Y' 



a 



n 



Kdr 



K 



+ -T^(n' + 2(x^-i)^ 



Kdr 



K 



1 

a 



1 2 



a. 



{2-^)iX'-lf 



Kdr 



(107) 



Integrating the F'-term and using the boundary conditions 



fx = nnrj^ + iirj^ 



a 



X' 



ne^XY 
K 



1 

+ - 



nY' a ^ 



for Y gives 

1 2 

-r 



+ 2{l--){X'-iy}Kdr 



:io8) 



Note that for a less than the critical value of 8 the mass per unit length is greater than 
Tmrj"^ whereas for a greater than 8 the final term in the integral is negative so that one 
cannot immediately tell if /i is greater than or less than irnr]'^. 

By adding equations (|101|) and (|102|) , then multiplying by K, and using equation (^) 
to integrate by parts we find that 

1 



271 / {pr+p^)e'^Kdr 



4 



(A' 



'Kdr. 



(109) 



This will be a useful equation in the next section when we relate the angular deficit 
to the physical properties of the string. It is also worth remarking at this point if we 
integrate Pr +P4, over the four dimensional region U bounded by the hyperplanes t = to, 

t = to + 1 and z = zq, z = Zq + 1 we find 



/ {Pr + P^) v^F^d^x = 2n [ {pr + p^)e^Kdr 
Ju Jo 



2tt 
0. 



°° d 
dr 



dr J 



(110) 
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This result explains why in the flat space limit the integrated radial and tangential 
pressures (a and b) are equal and opposite, and hence the constant c in equation (^) 
in fact vanishes in fiat space. It also explains why the quantity defined by ( |109|) is very 
small in the weak: field limit. This result is also consistent with the notion that in a 
certain sense Cosmic Strings have no active gravitational mass since the above result 
shows that the mass per unit length defined using Tolman's expression for the mass per 
unit length vanishes 



I {T',-Tl-Ti-T!)^)d' 

POO 

2tt / {pr + p^)e^Kdr since Pz = 
Jo 



= 0. (Ill) 

We now examine the thin string limit of the energy-momentum tensor. As in the 
previous section we first transform to {t,x,y,z) coordinates where x = rcos0 and 
y = r sin (p. Then because takes the form 

^e{x,y) = \a{x,y) (112) 

we can use the results of section 2 to deduce that is associated to a delta function 
and 

[a,{x,y)]^l26^'\x,y) (113) 
where /i is given by ( |106|) [or (|108|) ]. 



As in the case of a Cosmic String in Minkowski space the energy-momentum tensor is 
associated to one of the form ( [59| ) but because of ( |11CI| ) c ~ and d unless the 



integral of Pr (and hence p<^) vanishes. In general this does not occur except in the 
case of critical coupling (see below). Therefore although the energy-momentum tensor 
is always defined as an element of the Colombeau algebra Qs it is only associated to a 
distributional energy-momentum tensor for the special case of critical coupling. 

4-2. Critical coupling 

The field equations for a Cosmic String coupled to gravity simplify considerably in the 
case of critical coupling (which with our choice of numerical factors is given by a = 8). 
With this value of a the final term in ( |108| ) vanishes and we obtain the expression 



oo 

/x = W + 7rr/M ^^e-^'^X'-'-^^-^) +-[-^-^iX''-l)\ ^Kdr. (114) 



Since the integral is positive this gives a bound for the mass as 

fi ^ Tcnrj'^. (115) 
Clearly this bound will be saturated if 

A{r) = (116) 



ne^XF\2 1 



nY' a 
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and X{r) and F(r) satisfy the curved space Bogomol'nyi equations 
dX XY 

-d^ = ^^ ^'''^ 

^ = Ak{1-X^). (118) 
dr n 

It is readily verified that if in addition K{r) satisfies the equation 

= -8nr]^ \2K\X^ - if + n^X^Y^] (119) 

dr 

[i.e. equation (^) with A = and a = 8] then X(r), F(r), A{r) and K{r) satisfy 
the field equations for a critically coupled Cosmic String. If one imposes the boundary 
conditions at the origin 

X(0) = 0, r(0) = 1, K{<d) = 0, K'{0) = 1, (120) 

then one can look at the asymptotic behaviour of the solutions to ( 117 ) and (|118|) for 
large r and show that in addition 

lim X(r) = 1, lim Y{r) = 0. (121) 

r— »oo r— >oo 

So that X{r), Y{r), A{r) and K{r) given by the above equations satisfy both the field 
equations (^)- (l85|) and the boundary conditions (^)-(|89|) for a critical Cosmic String. 
Conversely one can show that a solution of the critically coupled Cosmic String field 
equations is also a solution of the curved space Bogomol'nyi equations with A{r) = 
and K{r) a solution to (|119|) (see Linet and Comtet and Gibbons |]T7|). 

It now follows from ( p.01| ) and (|102|) that the transverse pressures vanish so that 

Pr =0 (122) 

P^ = (123) 
and from ( |114| ) that 

/i = vrnr^l (124) 

It is quite remarkable that the mass per unit length for a critically coupled Cosmic String 
in a curved spacetime does not depend upon the details of the geometry but only on the 
winding number and t]. Because the transverse pressures vanish, the energy-momentum 
tensor of the 1-parameter family of Cosmic Strings (which gives the thin string limit) is 
associated to an ordinary distributional energy-momentum tensor. So that 

[T^^,{t,x,y,z)] ^W5(2)(x,i/)diag(l, 0,0,1). (125) 

Thus, as in Minkowski space, a critically coupled string can be modelled in the thin 
string limit without worrying about the internal structure (since it has a description in 
terms of classical distributions) but that for a non-critical Cosmic String the internal 
structure can not be completely neglected but can be described within the Colombeau 
algebra (which has a finer structure than classical distributions). 
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5. Angular deficit, iiolonomy and gravitational lensing 

For the conical metric given by equation (0) the spacetime is locally equivalent to 
Minkowski space as may be seen by changing to the new angular variable (p = ^0- 
In terms of this variable the metric becomes 

ds^ = dt^ - dr"^ - - dz^ . (126) 

However the spacetime is not globally fiat since 0^0^ 2ttA. Thus ( |126|) represents 
Minkowski space with a wedge of angle 27r(l — A) removed and the edges identified. 
In a similar way the Cosmic String spacetimes considered in the previous section 
are not asymptotically fiat, but are asymptotic to conical spacetimes. A convenient 
mathematical way of measuring the angular deficit is in terms of holonomy. If one 
parallel propagates an orthonormal frame once around the axis in a conical spacetimes 
one finds that it has been rotated about the axis by the deficit angle A0. Thus the 
holonomy gives an alternative way of measuring the angular deficit. For the Cosmic 
String spacetimes we will calculate the holonomy given by parallel transport around the 
curve r = ro, lying in the 2-surface t =const. z =const. and then look what happens as 
ro oo. This will give us the angular deficit of the limiting conical spacetime. 

Let the metric be given by 

ds^ = e^idt^ - dr^) - r^e^'dcj)^ - e^dz^ (127) 
and let be the orthonormal tetrad given by 

e^ = ?, e^ = P, e^ = 0^ = F (128) 

where the frame indices are in bold font to distinguish them from the coordinate indices. 
We now consider parallel propagation around the circle 

t = to, r = ro, 0e[O,27r], z = zq. (129) 

The unit tangent vector to this circle is given by 62 and if 7^^, denote the Ricci rotation 
coefficients then one may show that 

Ve,ei = 7?2eb = 7?2e2 = r-'e^^^^y^re^^'ye, (130) 
Ve,e2 = 72'2eb = 72'2ei = -r-'e^^^^y'{re^/'ye^. (131) 

So that a vector in the ei A 62 plane remains in that plane under parallel propagation 
around the circle. Let X = p(0)ei + g(0)e2 be a vector in the ei A 62 plane then from 
the above 

Ve,X = (e^p,^)ei + pVe,ei + (e^g,^)e2 + qV.^e^. (132) 



But 



and thus 



r 



(134) 
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Hence VeaX 



= implies 

dp 

d(j) ^ 
dq 

— = —Kp 



where 



-A/2 



d 

dr 



r=ro- 



(135) 
(136) 

(137) 



Solving this equation gives 

p{2n) = cos{27tk)p{0) - sin(27rK)g(0) (138) 

g(27r) = sin(27rK)j9(0) + cos(27rK)g(0) (139) 

so that the vector X rotates by — 27rK relative to the frame ei, 62- On the other hand 
this frame rotates by 27r as it goes round the loop, and X is rotated by 27r(l — k) as a 
result of parallel propagation around the circle r = tq. The holonomy at infinity gives 
the angular deficit and hence 



A0 

or in terms of A and K 



lim 2iT 

r—*oo 



lim 277 

r— +00 



(re 

dr 



d 



B/2 



dr ^ ' 



(140) 



(141) 



Because the t =const. z =const. 2-plane has zero extrinsic curvature the holonomy 
is also the same as that produced by the connection of the induced 2-metric. The 
holonomy in two dimensions may also be computed using the Gauss-Bonnet theorem 
(see eg Vickers [Q) and is given by 



Hol{Dr) 



2 — ^2 



(142) 



where Dr is the disk of radius r, fi^g is the curvature 2- form in the 2-plane and f2^2 
is the curvature 2-form of the 2-plane in the spacetime (which are the same since the 
extrinsic curvature vanishes). Taking the limit as r — > cxd gives 



A0 



2- 



Using the result of Proposition 3 in the Appendix we obtain 
A0 = Svr/i + 27r(P, + P^) 

where 



and 



Pr = 2tx I prKdr, 
'0 



= 27r / p^Kdr. 
10 



(143) 



(144) 



(145) 



(146) 



The thin string limit of Cosmic Strings 



20 



Note that equation ( |109| ) shows that this is consistent with Garfinkle's result (equation 
(82) of [^). Using equations ( |100| ), ( |101D and ( |102| ) we may obtain an exphcit expression 



for A0 in terms of the matter fields similar to equation ( |106|) for /i. 



A(j) = M'^rf 



AldXV 2 A^^y^ ,^2 (dY\ 



Kdr. (147) 



dr ) ' ' 2aK^ \dr J 

However the significant thing for us is that ( |144| ) gives an expression for the angular 
deficit in terms of the integrated matter fields. 

If one now considers the 1-parameter family of metrics given by Proposition 2, which 
give the thin string limit, one finds that all the terms in (|144|) are independent of e so 



that (|144|) is true for both the physical spacetime and the thin string limit spacetime. 
One may therefore give an expression for the angular deficit of the physical spacetime 
in terms of the distributional energy- momentum tensor of the thin string limit. Note 
that this result is not incompatible with Futamase and Garfinkle's claim that the thin 
string limit is not well defined using conventional distribution theory since the and 
terms are not well defined unless one uses a theory such as Colombeau algebras which 
extends the concept of distribution. 

The calculation above shows how the deficit angle A0 can be obtained in terms of the 
matter distribution. However the deficit angle can not be directly measured so we now 
consider the defiection of a light ray due to the gravitational field of a Cosmic String. 
One may calculate this by looking at the equations for the null geodesies. For the 
metric given by ( p7| ) these may be calculated using the Euler-Lagrange equations for 
the Lagrangian 

£= i(e'^f2-ey-pV02-e'^i2) (148) 

We will consider geodesies moving in the plane orthogonal to the string so that z = const, 
and z = 0. We now use the absence of explicit r and dependence to show that we 
have two conserved quantities: 
dC 

= e^T = E = const. (149) 

OT 

dC 

= -p^eV = J = const. (150) 

dcj) 

One may then calculate the scattering angle for a pair of geodesies which pass either 
side of the string. One finds that in these coordinates the angular defiection (pd is given 
by 

AE dp , , 

^^ = '""-1 MW)- - 11"^ 

where g{p) = j-e^^''~"'^ and pc is the closest approach of the light ray to the p = axis 
(which satisfies the condition g{pc) = !)• 
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In general we see that the scattering depends upon I = E/J which is essentially just 
the impact parameter for the geodesic. However for the special case of a conical metric 

ds"^ = dr'^ - dp^ - A^p^dcj? - dz^ (152) 

One can show that the scattering is independent of /. More precisely one finds that 

4 r ^dr 27r(l-A) 

If one now transforms to the coordinates one finds 

0rf = = 27r(l - A) = A0 (154) 

So that in the coordinates in which the cone is written as Minkowski space minus a 
wedge, the scattering angle is equal to the angular deficit. 

If one considers a possible observation of a multiple image of a star due to gravitational 
lensing by a Cosmic String then the closest approach of the geodesic would be 
approximately given by pc ~ -RA0, where R is the distance to the star. Thus the 
ratio of the closest approach to the thickness pd p would be very large. This implies 
that the impact parameter would be large compared to the thickness of the string. Thus 
when looking at the effect of gravitational lensing in the thin string limit the appropriate 
limit to consider is to take e — while keeping the impact parameter /, constant. We 
call this the lensing limit. In calculating the lensing limit the deflection takes place in 
a region which gets closer and closer to a conical spacetime as e — > 0. One can use this 
to show that in the lensing limit the scattering is independent of the impact parameter 
and is equal to that of a cone with deficit angle A0 given by equation (|144|) . 

However as well as the lensing limit one can also examine the limit in which the closest 
approach pc scales with the thickness of the string so that pc/ p is independent of e. 
We call this limit the grazing limit. In this case the scattering angle turns out to be 



independent of e and is given by ( [L51| ) for all values of e. The deflection in this case does 
of course depend upon the impact parameter, but a suitable choice to consider is / = 1 
since this gives the angle between to geodesies which are initially parallel and with a 
separation approximately equal to the effective diameter of the string. The results are 
plotted in Figure 1 which show 0^ and A0 as a function of a and rj respectively. 
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Figure 1. Scattering angle compared to deficit angle as a function of a and 77 
respectively. Notice that <j)d — A0 for a = 8. 



6. The Conical approximation 



We now wish to consider the asymptotic behaviour as e ^ of a 1-parameter family of 
solutions representing the thin string limit. As a result of proposition 2 we know that 
both the string and metric variables scale according to 

Mp) = f{p/e). (155) 

Thus in the 'far zone' given by p > e the asymptotic behaviour of /e(p) is given by the 
asymptotic behaviour of /(r) for large r. This gives the leading order behaviour of the 
metric variables as 

A{r) ~ Ao (156) 

B{r) ~ Bo (157) 

where Aq and Bq are constants such that 

ds^ = (dt^ _ dr^) - r^e^'dr^ - e^°dz^ (158) 

is the metric of a conical spacetime with deficit angle A0 given by (|144|) . Calculating 
the higher order asymptotic behaviour shows that the metric variables have asymptotic 
expansions of the form 

A{r)=Ao + ^ + 0{^) (159) 

5(r) = 5o + ^ + 0(i). (160) 

In contrast to the asymptotic expansion of the metric variables in inverse powers of 
r one finds that the matter variables have exponentially fast decay. If one defines 
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Z{r) = X{r) 
equations 



1 and W{r 
d 



Y{r)/r then to leading order X(r) and Y{r) satisfy the 



dr 
d 



dr J 
dW 



dr \ dr 



1 + ae^°r2)Vr = . 



(161) 
(162) 



These are hyperbohc Bessel equations of order and 1 respectively and give the 
asymptotic behaviour as 

Z{r) ~ Ko{cir) ~ r^^/'^e-^^' (163) 
W{r) ~ Ki{c2r) ~ r-^/^e'^^^ (164) 

where ci = 2^/2e^°/'^ and C2 = ^/ae^°/'^. 

Calculating higher order terms leads to asymptotic expansions for X[r) and Y{r) of the 
form 



Xir) 



Yir) 



l + r-i/2g-cir 



Xo + ^ + 0(l) 



lo + - + 0(4) 



(165) 
(166) 



Because of the rapid exponential decay of the matter terms compared with the 
polynomial in 1/r decay of the metric variables, one finds that even in the 'near zone' 
(p < e) provided one remains outside the exponentially small inner 'core zone' then 
the asymptotic behaviour of the metric is the same as that of a cylindrically symmetric 
vacuum spacetime which is asymptotic to a conical spacetime. We therefore introduce 
the 'conical approximation' in which we model the matter as a conical singularity with 
angular deficit given by ( |144| ) and consider the gravitational field as a perturbation of this 
metric. Because the gravitational field satisfies a linear field equation in this situation 
one can solve for the perturbation exactly in terms of a cylindrically symmetric solution 
of Laplace's equation and then write the gravitational field as a superposition of this 
with the conical solution. 



The asymptotic behaviour for the above variables for large r translates into the following 
thin string limit outside the core. 

(167) 



Up) 

Se{p) 
Peip) 



P 



&0 + — + 0(e^) 
P 



V 



^l/2^-l/2g-c,pA 



^-l/2^1/2g-c,pA 



So + — + O(e^) 
P 

P 



(168) 
(169) 
(170) 
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7. The thin string hmit for dynamic Cosmic Strings 

In this section we consider the thin string limit for a time dependent Cosmic String. We 
start with the Lagrangian (|^) and with $ and given by (|T^) and (|TT]) but now take S 
and P to be functions of both r and p. We first consider the case of a dynamic Cosmic 
String on a fixed Minkowski background. The equations of motion are then given by 

We now introduce the Colombeau thin string parameter e and use this to scale the 
coupling constants e, A and rj according to equations (|29|)-(|3TD. We then consider 1- 
parameter families of solutions and to the equations 

% (^^) - = '^^''^'^^'^ - + ^'''^ 

The solutions to ( |17ij| ) and ( |174D are given by 

S,{r, p) = r7X(v^77p/e, ^/Xr]r/e) (175) 
P,(r, p) = Y{VX7]p/e, VX7]T/e) (176) 

where X{t,r) and Y{t,r) are solutions to 

4(r^) - '-'S = •^|4r^(A- - 1) + (177) 



9r V dr ) dt^ 



9r V dr ) dt"^ 

We see from the form of P^ and given by ( |175| ) and ( |176| ) that taking the limit as 
e — > corresponds to simultaneously taking a limit in which both p oo and r ^ oo. 
Since a dynamic Cosmic String eventually settles down to a static configuration taking 
the thin string limit is also a limit in which the string also becomes static. 

In order to investigate the asymptotic behaviour for small epsilon (for p 7^ 0) we consider 
solutions of ( |1 771 ) ^"^^ (|1 78|) for large values of r. These may be written 



X = 1 + Xi (179) 

Y =Yi (180) 

where Xi and Yi are small for large r. Substituting into ( |177|) and ( |178| ) and ignoring 
terms which are quadratic or higher in Xi and Yi we find they satisfy 
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We may solve ( |181| ) by considering a solution of the form 

X^{t,r) = Z{r)e"^' (183) 
and substituting into ( |181|) to obtain 



f) / r)7\ 

r^ir^] -i^X-A^Z = Q . (184) 
dr\ or J 

The solution to this which is finite as r ^ oo is given by Z{r) = Ko^cir) where 
c? = 8 — cu^. Hence 



X^{t, r) = Ko(V8^TJ2r)e'"* (185) 
and the general solution to ( |181|) may be written 



X{t,r) = J A^KoiVS - u^r)e''^'du . (186) 

In terms of the physical variables this gives the asymptotic behaviour in the far zone 
(p > e) as e — > as 

5,(p, r) ~ 77 + y A^iro( V(8A - uj^)X7]^p/e)e'^^^'^^/^du; . (187) 

In a similar way one can show that 

Pe{p,r) r-^P-j B^K^{^{eW - uj^)W p/ t)e'^^''^^^l'duj . (188) 

The important thing to note is that 5*^ and still have an exponential decay as e ^ 
due to the exponential decay of Kq and Ki for large r so these expansions remain valid 
in the near zone provided one stays outside the exponentially small inner core. 

A similar analysis of a Cosmic String coupled to the gravitational field shows that and 
Pe again have the e~'^pl'^ type exponential fall off. By contrast the metric coefficients have 
asymptotic expansions which are polynomial in e. In the dynamic case we may therefore 
also introduce the conical approximation in which we model the matter as a static conical 
singularity with angular deficit given by (|144|) and consider the gravitational field as a 
perturbation of this solution. Because of the linearity of the vacuum Einstein equations 
for a cylindrically symmetric metric written in Rosen form, one can explicitly write down 
the solution as a superposition of a conical solution with a solution of the cylindrical 



wave equation. Such a situation has been analysed by Marder |19|. 



Another important feature of (|187|) and ( p.88| ) is that they may be used to determine 



the frequency of oscillation of the matter fields. Since K^) and K\ decay exponentially 
fast, the dominant contribution to the frequency is when the argument is close to zero. 
This gives the frequencies of the fields S and P as 

fs = V2Xr]/7T (189) 
fp = er]/27T . (190) 

A plot of the corresponding values for the rescaled fields X and Y, compared to the 



numerically computed values (using the results of |^) is shown in Figure 2 
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Figure 2. The frequecies fx and fy calculated using the asymptotics compared to 
the numerically computed values. 



8. Conclusion 



In this paper we have examined the thin string hmit of a Cosmic String coupled to the 
gravitational field. By using a description in terms of the simplified Colombeau algebra 
we have been able to show that the energy- momentum tensor has a well defined limit. 
Furthermore we have shown how the angular deficit of the limiting conical spacetime 
may be expressed in terms of the limiting energy-momentum tensor through equation 
( |144|) . However it is important to note that in general the Colombeau distributional 
energy-momentum tensor that one obtains is in general not associated to a conventional 
distribution. It is only in the case of critical coupling that the energy-momentum tensor 
can be written in terms of Dirac delta functions. Thus the results we obtain are not in 
conflict with those of Garfinkle for example 0] who showed that there were problems 
describing the thin string limit. It is only by using the features of the Colombeau algebra 
which allow one to describe the internal structure of the line source that it is possible 
to describe the thin string limit. For example a limiting distributional tangential stress 
cannot be defined using conventional distributions but can be described in terms of the 
Colombeau algebra. 

In the last two sections we examined the asymptotic behaviour of both the metric and 
matter variables in the thin string limit. We showed that in the far zone the metric 
variables have polynomial decay while the matter variables decay exponentially. We 
were therefore able to introduce in both the static and dynamic case the notion of the 
conical approximation. In this approximation one replaces the matter terms by a line 
source with distributional energy-momentum tensor and a conical singularity which has 
angular deficit given in terms of the matter through (|144|) . One is then able to model 
the gravitational field as a solution of the vacuum Einstein equations with the required 
angular deficit on the axis. Because of the linearity of the field equations one can 
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solve this explicitly by looking at cylindrically symmetric solutions of the wave equation 
which are regular on the axis and looking at a superposition of this with a static solution 
representing a conical spacetime. Because of the exponential fall off of the matter fields 
this approximation is valid even in the near zone as long as one remains outside the 
exponentially small inner core of the string. We have shown how this approximation 
may be used to calculate the frequency of oscillation of the matter fields and this shows 
very good agreement with the numerically calculated frequencies of |20|- ^ future 



paper this idea will be developed and the conical approximation will be used to examine 
further properties of a dynamic Cosmic String and the results compared to numerical 
calculations. 

9. Appendix: A formula for the deficit angle in terms of integral quantities 

In this appendix we obtain an expression for the angular deficit in terms of integrals of 
the components of the energy-momentum tensor. The calculation just uses the Einstein 
tensor and is entirely in terms of the metric coefficients. No use is made of the field 
equations for the matter variables. In this section we take the metric to be given by 

ds^ = - e^'dr^ - e^'d4>' - e^'^dz' (191) 

(Note the factors of 2 in the exponents compared to earlier sections.) and let e^ be the 
orthonormal tetrad given by 

e^ = ?, e^ = P, e^ = ?^ e^ = F (192) 

where the frame indices are in bold font to distinguish them from the coordinate indices. 
Then the calculations of section 5 show that the angular deficit is given by 

A0= lim 27r(l-cV"^). (193) 

r— +00 

Using the Gauss-Bonnet theorem one can also write this as an integral 

/■oo 

A0 = -2tx / [c" + {df - c'6']e^~Mr. (194) 



^0 

On the other hand using a = Too = l/(87r)Goo (where the use of bold font indicates the 
use of frame indices) we get 

a = -l/(87r)[a" + c" + {a'f + (c')' - b'd - a'b' + a!d]e-^\ (195) 

So that the mass per unit length is given by 

r"27r fOO 



av^drd0 (196) 

^0 
/■oo 

2tx I (xe^+'^dr (197) 

+ [df - b'd - alb' + a'c'le^-Mr. (198) 







1 ["^ 

-J [a" + d' + {a 
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poo poo 

/ a"e^-Mr = [a'e^-^]^ - / {a'c' - a'ty-'^dr (199) 
Jo Jo 

poo 

= - / (a'c' - a'ty-'^dr. (200) 
Jo 



Hence 



and thus 



1 7"°° 

pL^-^ [c" + {c' f - b'c' + (a')2]e^-Mr (201) 



4^0 



pco poo 

8nn = -27r / [c" + (c')^ - 6'c']e"-Mr - 27r / (a')2e'^-Mr (202) 
io Jo 

poo 

= A0 - 27r / (a')2e"-Mr. (203) 
Jo 

We now show how to write this entirely in terms of integrated components of the energy- 
momentum tensor. Using — Tii = l/(87r)G'ii we have 

Pr = l/(87r)[(a')2 + 2a'c']e-2'' (204) 

Let Pr be the integrated version (with respect to the volume form determined by the 
determinant of the 2-metric on t=const. z=const.) Then we have 

p2n poo 

Pr^ / Pre'+Mrd0 (205) 
Jo Jo 

= ^ y [(a')2 + 2a'c']e^-Mr. (206) 

Similarly using = T22 = 1/ (871)^22 we have 

p^ = l/(87r)[2a" + 8(0')^ - 2a'b']e-^^ (207) 

so that the integrated version is given by 
1 f°° 

P<t>^- [2a" + 3(a')2 - 2a'6']e^-Mr. (208) 
4 Jo 

But (as before) 

/*CXD poo 

/ a"e^-Mr = [a'e^-'']^ - / (a'c' - a'6')e'^-^dr (209) 
Jo Jo 

poo 

= - / (a'c' - a'6')e"-Mr (210) 



so that 



Thus 



So that 



1 /""^ 

P^^- [3(a')2 - 2a'c']e'^-Mr. (211) 
4 Jo 

/■oo 

- / {a'fe^-^dr. (212) 
io 



000 

Pr + -P^ 



A(/) = 87r^ + 27r(P^ + P^) . (213) 
We have thus established the following proposition. 
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Proposition 3 The angular deficit A(f) for the metric given by 
ds" = e^'^df - eV - e^^d^^ - e'^'dz^ 

where a, b and c are only functions of r is given by 
A(j) = Swfx + 27r(P, + P^) 



(214) 



(215) 



where ji, Pr and P^ are the integrated values o/ 1/ (871)^00; l/(87r)G'ii and 1/(871)^22 
respectively, and the integration is with respect to the volume form induced on the 2- 
surface t = const, z = const. 
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